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Abstract 


In 2016, an exposure meter was installed on the Lijiang Fiber-fed High-Resolution Spectrograph to monitor the 
coupling of starlight to the science fiber during observations. Based on it, we investigated a method to estimate the 
exposure flux of the CCD in real time by using the counts of the photomultiplier tubes (PMT) of the exposure 
meter, and developed a piece of software to optimize the control of the exposure time. First, by using flat-field lamp 
observations, we determined that there is a linear and proportional relationship between the total counts of the PMT 
and the exposure flux of the CCD. Second, using historical observations of different spectral types, the 
corresponding relational conversion factors were determined and obtained separately. Third, the method was 
validated using actual observation data, which showed that all values of the coefficient of determination were 
greater than 0.92. Finally, software was developed to display the counts of the PMT and the estimated exposure 
flux of the CCD in real-time during the observation, providing a visual reference for optimizing the exposure time 


control. 
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1. Introduction 


The signal-to-noise ratio (S/N) is one of the key parameters 
to evaluate the data quality of high-resolution spectra. During 
observations, because we cannot obtain the CCD exposure flux 
information directly, the exposure time is generally estimated 
by considering the weather and observing conditions, and the 
characteristics of stars, etc, which have large uncertainties. If 
the exposure time is blindly increased for high S/N, it is a 
waste of observation time, and it may pollute the observed data 
because the CCD is saturated. If the S/N does not meet the 
requirements, repeated observations will affect the observation 
efficiency. There are usually two methods to solve the problem. 
One is using the Exposure Time Calculator (ETC), which needs 
to accurately simulate the target, weather and instrument 
information before observations. Most instruments use ETC to 
estimate exposure time, although it has some uncertainties. The 
simulated S/N from the ETC is overestimated by 40%-50% on 
IGRINS (Le et al. 2015), while it is within 10% for HARPS 
(ESO 2021). 

Another way is using Exposure Meter (EM), which is 
installed as a part of the instrument. the EM can monitor the 
flux in real-time to optimize the exposure time, because it can 
terminate the exposure when the desired counts have been 
acquired. The EM can also calculate the photon-weighted 
mean time of each exposure, which is critical for precise 
RV(Radial Velocity) measurements (Tokovinin et al. 2013; 


Telting et al. 2014). For precise RV measurements with long 
exposures, EM even becomes a necessity (Raskin 2011). Many 
high-resolution spectrographs have been equipped with EM, and 
the detectors of EM are photomultiplier tubes (PMT) or CCD 
(Kibrick et al. 2006; de Cuyper et al. 2007; Gibson 2013; 
Tokovinin et al. 2013; Landoni et al. 2014; Ben-Ami et al. 2016; 
Blackman et al. 2020; Gupta et al. 2022). EM usually picks up 
from 1% to 8% of the light from behind the slit or grating. 
According to the different spectrographs, the count rates of EM 
are different. For a star of V=5 mag., the typical count rate is 
103 s_' in the EM of CHIRON (Tokovinin et al. 2013), although 
its maximum count rate is 456 ae and the error of the mean 
exposure time is within 1s. However, for a star of V = 16 mag., 
the count rate is 19 s~', and the error of the mean exposure time 
is within 165 s for the EM of HARPS (Chile 2010). 

In this paper, we propose a method to estimate the exposure 
flux of a CCD in real-time using the counts of PMT in Fiber- 
fed High-Resolution Spectrograph (HiRES) of the Lijiang 
2.4m telescope. Through flat-field lamp observation, we 
determined that the signals between the two detectors are 
linearly proportional, and obtained the relationship conversion 
coefficients for six different spectral types using historical 
observed data. The validation results from the observed data 
show that the estimated values of the CCD exposure flux are in 
good agreement with the actual observed values. At the same 
time, the application software was developed to realize the 
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Figure 1. The optical system of HiRES for 2.4 m Telescope, the red part is photonic detection system. 


real-time display of the estimated CCD exposure during 
observation. Section 2 describes the principle of the method. 
Section 3 shows the validation results, Section 4 describes the 
working principle of the software and shows the user interface, 
and Section 5 summarizes the results of this work and discusses 
the error factors. 


2. Linear Modeling for Exposure Meter 
2.1. Optics of HiRES and Exposure Meter 


HiRES of the Lijiang 2.4 m telescope has two science fibers, 
the spectral resolutions (R = A/AA ) are 32,000 and 49,000 for 
the 20 and 1”2 diameter fibers (on the sky) at 550nm, 
respectively. The orders of an observed spectrum are from 61 
to 156, and the wavelength coverage is from 380 to 990 nm. 
The long-term stabilities of temperature and pressure for the 
spectrograph are controlled within 26 + 0.25°C and 30 + | Pa, 
respectively (Wang et al. 2019). A starlight coupling device 
was installed on a side port in the AG-box unit (Fan et al. 2015) 
at the Cassegrain focus of the telescope and guided the light 
from the telescope or calibration lamps to the spectrometer via 
scientific fiber. Figure 1 shows the optical layout of HiRES of 
2.4 m telescope. 

An EM (red box in Figure 1) was installed in the 
spectrograph optical path for monitoring starlight coupling 
into science fiber during observation to pick up 3% starlight by 
a small folding mirror to EM, and its main component is a wide 
dynamic range PMT. Figure 2 shows the quantum efficiency 


curves of the two detectors (CCD and PMT), they have high 
detection efficiency in the wavelength range from 500 nm to 
800 nm, and the peak quantum efficiency of PMT and CCD 
detectors are 400 nm and 550 nm, respectively. We used the 
observed data to study the relationship between the counts on 
EM and the flux on CCD of spectrograph, and developed a 
piece of software to improve the observation efficiency. 


2.2. Mathematical Model 


We first investigated the relationship between the signals 
from the two detectors (CCD and PMT). This was carried out 
by using flat-field lamp observations. Specifically, different 
exposure times (0.5 s, 0.6 s, 0.7 s, 0.8 s due to the bright enough 
flat-field lamp) were set for group observations, and four sets of 
spectral images of the CCD and the corresponding total counts 
of the PMT were obtained. First, the background signals of the 
spectral images and PMT counts are corrected. Second, the 
2000th column of pixels with the strongest signal in the 
spectral image was taken, and the ADU values of the pixels 
covered by the 96 orders in this column were summed 
individually to obtain the exposure fluxes at the wavelength 
points corresponding to the different orders. Figure 3 shows the 
relationship between the exposure fluxes of the 10 orders and 
the corresponding total PMT counts. The horizontal coordi- 
nates represent the total PMT counts and the vertical 
coordinates represent the exposure fluxes of the different 
wavelength points in the column of the CCD image. 
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Figure 2. Distribution of detection efficiency of PMT (model: Hamamatsu H8259-02) and CCD (model: E2V CCD203-82). 
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Figure 3. The linear fitting between the observed flux on CCD and the counts of PMT, the scatters represent the sampled values of the observed flux of different 


exposure time. 


The dark counts for the CCD and the PMT were obtained 
before observations. In the no-exposure state, the effective 
signals of both detectors are zero. The relationship can be 
expressed in the form of Equation (1). Here Fecp is the 
observed flux at the corresponding wavelength point for each 
pixel of the CCD, Cpmr is the sum of counts of the PMT during 
the exposure time T, G is the scale factor matrix of the 
corresponding pixel on CCD image, dy is the dark of the CCD 
during the exposure time T, and by is the sum of the dark 


counts of the PMT during the exposure time T 
Fecp — dr = G(Cpmr — br). (1) 


The total count of the PMT is the sum of the exposure fluxes 
of all the bands it detects. The energy distributions of stars of 
different spectral types are also different, and with the 
inconsistent response of HiRES in different bands, as well as 
the actual atmospheric absorption, it is difficult to obtain more 
accurate spectral response curves using the standard blackbody 
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Figure 4. Two-dimensional images of a flat and a star, respectively. The left is flat, while the right is a star (HD189733). 
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Figure 5. The distribution of ADU (normalized) in different wavelengths observed by stars of different spectral types on CCD (2000th column). 


radiation formula. The actual observed stellar spectra are 
different from those of the flat field lamps as shown in Figure 4. 
In addition, Figure 5 shows the actual observed energy 
distributions for stars of six spectral types, where the observed 
energies have been normalized. There are many absorption 
lines in the star’s spectrum, and the right image shows that the 
darker parts of the discontinuity are absorption lines. Based on 
these spectral characteristics of HiRES, we need to determine 
the respective conversion factor G matrices for the stars of 
different spectral types. 

With Equation (1), we determined and obtained six types (Type 
B, A, F, G, K, and M) of conversion factors, using historical 
observations. These observations include: CCD spectral image, 
PMT total counts, CCD dark, and PMT dark counts. Figure 6 
shows the matrix of conversion coefficient for a star of spectral 
type A. With the above method, we can estimate the exposure flux 
of the target spectrum using the real-time count of the PMT and 
obtain a simulated CCD spectral image. 


3. Validation 


We have validated the above method by observations from 
2022 July 16 to 2022 July 20. The information of the observed 
targets is given in Table 1, and the CCD dark and PMT dark 
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Figure 6. The conversion coefficient matrix for spectral type A. 


counts were obtained from the observations on the same day. The 
conversion factor (G) corresponding to the spectral type of the 
target to be validated was chosen, and the flux (Feca) of the CCD 
was estimated by Equation (1) using the counts of PMT obtained 
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Figure 7. Comparison of estimated and observed values of flux for stars with different spectral subtypes. The blue curves are the observations, and the orange curves 


are estimation. 


Table 1 
Information on Observed Targets for Validation 
Spectral Type of Stars for Validation Magnitude R 
B9.5V 5.665 mag. 0.964 
AIV 3.7 mag. 0.970 
F8 7.58 mag. 0.964 
G1IIV 7.457 mag. 0.967 
K1.501 —0.05 mag. 0.923 
M1.0V 10.03 mag. 0.990 


from the observation. The results were evaluated using a 
coefficient of determination between the estimated value of the 
CCD exposure flux and the actual observed value. The coefficients 
of determination were obtained using Equations (2), (3), and (4) 


R2 = E SSestimated (2) 
SS total 
SSTotal = S20; = Ymean) (3) 
i=1 
SSestimated = “0; = ne (4) 


i=l 


where R° is the coefficient of determination, y; is the observed 
flux of the CCD, Yestimatea 18 the estimated flux of the CCD, and 
Ymean 18 the average of the observed fluxes. If the estimate and 
the observation are in perfect agreement, the value of R? will be 
1. The third column of Table 1 shows that all R? values are 
greater than 0.92, indicating a high degree of agreement. 
Figure 7 shows a plot between the estimated and observed 
values of the exposure fluxes on pixels of the CCD, where the 
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Figure 8. Flow chart of observation software. 


blue curve represents the observed values and the orange curve 
represents the estimated values. The absorption line and the region 
between the two orders in the estimated values are ignored for 
clarity of comparison. By comparing the results of the above 
figure and table, it can be concluded that our method can estimate 
the exposure flux of the CCD using the real-time counts of the 
PMT in the observation, and the estimation results are reliable. 


4. Software 


We developed a piece of software that uses the EM to 
estimate the exposure time based on Python and the Qt 
Designer IDE of the Windows 10 system. 
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Options: 
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Figure 9. Observation software operation interface. 


Figure 8 shows the software workflow. The working 
mechanism is mainly performed through the communication 
of three threads. These threads are the main thread, the data 
acquisition thread and the spectra of the thread. The main 
thread receives the operations and sends the command to the 
PMT counter through the serial port. The PMT counter does as 
commanded and sends data back to the software. The data 
acquisition thread is responsible for receiving counts, verifying 
the data, and subtracting the background counts. The main 
thread receives the photon counts, computes statistics (includ- 
ing sample value, mean value, total value, and number of 
samples), and dynamically updates the statistical results. The 
spectra of the thread is responsible for estimating the current 
exposure flux at a given frequency (0.2 Hz) and updates in 
real time. 

Figure 9 shows the user interface of the software. Observers 
can select or fill in the star’s information in the operation 
interface according to the scientific requirements before 
observation. The PMT counts and the estimated exposure flux 
can be displayed in real-time during observation, providing a 
visual and reliable reference for the observation. 


5. Conclusions and Discussion 


We have presented a method for estimating the exposure flux 
of a CCD in real-time using the counts of the PMT of an 
exposure meter, and show the validation of the method using 
observational data. The validation results show that the 
estimates obtained using the method are very much in line 
with the actual observations. 

A linear and proportional relationship between the total 
counts of the PMT and the fluxes of the pixels corresponding to 
different wavelength points of the CCD spectral image has 
been determined by flat-field lamp observations, and the 
equation of the transformation relationship is given. The six 
spectral types’ transformation coefficient templates were 
determined from historical observations. The validation of the 
observation data from 2022 July 16 to July 20 shows that the 
coefficient of determination values between the estimated and 
observed values of the CCD exposure fluxes are all greater than 
0.92. It is shown that the method can be used to estimate the 
CCD exposure flux using real-time counts of the PMT in the 
observation, and the results are reliable. In addition, the 
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working mechanism and user interface of the method’s 
application software are shown. The software displays in 
real-time the PMT counts and the estimated value of the 
exposure flux at the CCD at that time, providing a visual 
reference for optimizing exposure time control. 

In Section 2, we determine a template of model transforma- 
tion coefficients for different spectral types using historical 
observations, and the key lies in the choice of data. First, the 
data were chosen to cover as many subtypes of the spectral type 
as possible. Due to precious observing time, we could not 
observe targets of all subtypes. Second, higher quality 
observations were selected, such as those with high S/Ns 
under stable observing conditions. The distribution of PMT 
observations can be used to determine whether the observing 
conditions were stable at that time and whether they were 
affected by cloudiness. Considering the above factors, the 
actual amount of selected observations is small. Third, there are 
different magnitudes for stars of the same spectral type. We 
found that the response distribution of spectral imaging on the 
CCD varies between targets with large magnitude differences, 
especially in the blue-end wavelength region, due to the 
atmospheric window and spectrometer efficiency. Therefore, 
we must also perform separate model transformation coefficient 
determinations for particularly bright standard stars of the same 
spectral type, or for particularly faint special targets. 

In addition, the effects caused by the sky-light background 
and the atmospheric absorption line are not considered in this 
paper. According to the long-term monitoring statistics of 
Lijiang Observatory, the sky-light background is between 16.5 
and 21 mag. arcsec * (Xin et al. 2020). This may introduce 
some uncertainties in the estimation of faint stars observed with 
long exposures. This factor must be taken into account if high 
accuracy is required. For example, we can make an observation 
of the sky-light background in the region of the target as a 
background correction before the observation. To a certain 
extent, it can reduce the influence of the sky-light background, 
but it also loses additional observation time and affects the 
observation efficiency. In addition, the intensity of the 
atmospheric absorption line of Lijiang Observatory varies with 
the seasons (Lu et al. 2021), and the flux at the wavelength it 
covers will also follow the influence, and the observer can 
avoid referring to the pixel area covered by this wavelength. 

Future plans include increasing the sample size of observa- 
tions to optimize model transformation coefficients to improve 
estimation accuracy, and improving software functionality 
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based on observation requirements. This research provides 
valuable insights for the application of the method on other 
spectrometers of the same type. It should be noted that the 
current implementation of the method is limited to bands with 
wavelength coverage less than 10000 A in HiRES for the 2.4 m 
telescope. However, with the selection of an appropriate 
detector, this method can be extended to other wavelength 
spectrometers. 
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